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Abstract. We study the dynamical behaviour of the interacting holographic dark 
energy model whose interaction term is Q — SH{XdPd + AcPc), where pd and pc are 
the energy density of dark energy and CDM respectively. To satisfy the observational 
constraints from SNIa, CMB shift parameter and BAO measurement, if Ac = Xd or 
Ad, Ac > 0, the cosmic evolution will only reach the attractor in the future and the ratio 
Pc/ Pd cannot be slowly varying at present. Since the cosmic attractor can be reached in 
the future even when the present values of the cosmological parameters do not satisfy 
the observational constraints, the coincidence problem is not really alleviated in this 
case. However, if Ac 7^ A^ and they are allowed to be negative, the ratio pd pd can 
be slowly varying at present and the cosmic attractor can be reached near the present 
epoch. Hence, the alleviation of the coincidence problem is attainable in this case. The 
alleviation of coincidence problem in this case is still attainable when confronting this 
model to SDSS data. 

Keywords: Dark Energy. 



The coincidence problem and interacting holographic dark energy 



2 



1. Introduction 

Observations suggest tliat the expansion of the universe is accelerating [1], [2]. The 
acceleration of the universe may be explained by supposing that the present universe is 
dominated by a mysterious form of energy whose pressure is negative, known as dark 
energy. One problem of the dark energy model is the coincidence problem, which is 
the problem why the dark energy density and matter density are of the same order of 
magnitude in the present epoch although they differently evolve during the expansion 
of the universe. A possible way to alleviate the coincidence problem is to suppose that 
there is an interaction between matter and dark energy. The cosmic coincidence can 
then be alleviated by appropriate choice of the form of the interaction between matter 
and dark energy leading to a nearly constant ratio r = pc/ Pd during the present epoch 
[3], m E] or giving rise to attractor of the cosmic evolution at late time P, [7] . Since the 
existence of the cosmic attractor implies constant r but the attractor does not always 
occur at the present epoch, we first find a range of dark energy parameters for which 
the attractor occurs and then check the evolution of r during the present epoch. 

Based on holographic ideas [H [9], one can determine the dark energy density in 
terms of the horizon radius of the universe. This type of dark energy is holographic 
dark energy [10] - [H]. By choosing Hubble radius as the cosmological horizon, the 
present amount of dark energy density agrees with observations. Nevertheless, dark 
energy evolves like matter at present, so it cannot lead to an accelerated expansion. 
However, if the particle horizon is chosen to be the cosmological horizon, the equation 
of state parameter of dark energy can become negative but not negative enough to drive 
an accelerating universe. The situation is better when one uses the event horizon as 
the cosmological horizon. In this case, dark energy can drive the present accelerated 
expansion, and the coincidence problem can be resolved by assuming an appropriate 
number of e-foldings of inflation. Roughly speaking, the coincidence problem can be 
resolved because the size of the cosmological horizon during the present epoch depends 
on the amount of e-folds of inflation, and the amount of holographic dark energy depends 
on the horizon size. Nevertheless, the second law of thermodynamics will be violated 
if < — 1 [H [15]. Hence, Wd should not cross the boundary Wd = —1- The boundary 
= — 1 can be crossed if dark energy interacts with matter. Since now the horizon size 
has a dependence on the interaction terms,the alleviation of cosmic coincidence should 
also depend on the interaction term. 

In this work we suppose that the holographic dark energy interact only with cold 
dark matter (CDM) and treat baryons as non-interacting matter component. Our 
objective is to compare the region of dark energy parameters for which the cosmic 
evolution has an attractor within the parameter region that satifies the observational 
constraints from combined analysis of SNIa data [I6], CMB shift parameter [17] and 
BAO measurement p^. The results of the comparison can tell us about the range of 
parameters that alleviate the cosmic coincidence. 
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2. The autonomous equations 

In this section, we derive the first order differential equations that describe the evolution 
of radiation, baryon, CDM and dark energy densities in the universe. By analyzing these 
equations, one can estimate the asymptotic evolution of the universe. To proceed, we 
start from the Friedmann equation 

H^ + ^ = -^{Pr+Pb + Pc + Pd) , (1) 

d oTTlp 

where H is the Hubble parameter and the subscripts r, b, c and d correspond to the 
radiation, baryons, CDM and dark energy respectively. The parameter K denotes the 
curvature of the universe, where K = — 1,0, +1 for the close, fiat and open universe 
respectively. The above equation can be written in terms of the density parameters 
Qk = K/{a'^H^) and fi„ = pa/iSrhlH^) as 

l + nK= = fir + fib + fic + fid- (2) 

a=r,b,c,d 

The index a runs over the 4 species, namely radiation, baryons, CDM and dark energy. 
We now derive the autonomous equations for the dynamical variables Qk and fia. 
Differentiating fi^ = K/i^a^H"^) with respect to In a, we get 

where prime and dot denote derivative with respect to In a and time respectively. From 
the definition of the density parameter, one can show that 

^'. = ^A^-^^\- (4) 




Hp, 

To study the evolution of the universe at late time, we will search for the fixed points 
of the above autonomous equations and check the stability of these fixed points. The 
fixed points of eqs. ([3]) and (j4]) are the points (fii^c, fiac) at which 

^'K = ^'a = 0. (5) 

It follows from eq. ([3]) that fi^ = at VL^ = or 1 + H/H^ = 0. Hence, possible fixed 
points at Qk 7^ correspond to the non-accelerating universe, i.e. 1 + H / H"^ oc a = 0. 
Since the expansion of the universe is accelerating today, we consider only the fixed 
points at fix = 0, and therefore neglect fix in our consideration for simplicity. 

Now we come to the case of interacting holographic dark energy and use eq. (jl]) to 
obtain the autonomous equation for this case. In the holographic dark energy scenario, 
the energy density of dark energy is related to the cosmological horizon L by 

Pd = Sc^mpL"^ (6) 

where c is a positive constant. Differentiating the above equation with respect to time, 
we obtain 

Pd = -^Pd^- (7) 
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We take the cosmological horizon to be the event horizon, which is defined as Re{t) = 
a(t) J^di/a(t). Hence, 

Li rie -iXe 

We therefore get 

3/2 

p, = -2Hp, + 2^^. (9) 
When dark energy has an interaction with CDM, the continuity equations yield 



Pc = -3Hp, + Q, (10) 

Pd= -3H{1 + Wd)pd-Q, (11) 

where Q = SH^XdPd + KPc)- Usually, one supposes that Q > because the second 
law of thermodynamics might be violated if energy transfers from matter to dark 
energy {Q < 0). However, for generality, we will not restrict Q to be positive in our 
consideration. Comparing eq. ([9]) with eq. ffTTl) . we get 

Wd = ---2— . (12) 

6 6C \Ld 

We assume that radiation and baryons have no interaction with dark energy, so that 
they obey the continuity equations 

Pr = —AHpr and ph = —3Hpb. (13) 

From eq. ([I]), one can show that 

2th = ^-FFT = -^ + ^d + 2^ + 3(A,^], + AA) - (14) 
In the above equation, we set Qk = 0. Using eqs. ([9]), f|T3|) and f|T^ . we obtain 



n'^ = nAi + 2^ - firf - 2^ - 3(Arfnrf + A A) + a, ) , (is) 

/ 3 fi^/^ 
^', = ^ci -^i>^d^d + -^d- 2^ - 3(AA + XM + a. 1 , (16) 

= a i-l-^d- 2-^ - SiXd^d + + a. 1 , (17) 

The fixed points of the above equations are {Qd, ^c, ^r, ^b) = (0, 0, Qrc, 0), (0, 0, 0, Qbc) 
and {fide, flee, 0, 0). Since we are interested in the late time evolution of the universe, we 
will consider only the fixed point {fide, flee, 0, 0). This fixed point can occur at late time, 
i.e., about the present or in the future, because flr/flc usually decreases with time and 
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fifo/i^c can decrease with time if Q > 0. Using eqs. (fTBI) and (ITB]) . the relation between 
Qdc and Qcc can be written as 

1 + 2^ = —{X,n,, + Aefiec). (19) 

From eq. ([2]) , we have fi^c + Qcc = 1 and hence 

1 - 3A, + - (1 + 3Xd - 3Xc)ndc - = 0. (20) 

c c 

The solution of eq. (l20l) gives fi^c in terms of A^, Ac and c. Instead of finding the solution 

of this equation, we will use eq. (|T9l) to compute the cosmological parameters of interest 

in the following section. The basic idea is that there are various values of c, Xa and Ac 

that satisfy eq. f|T9|) for given Q^c and Qcc- Changes in the values of c, Xd and Ac lead 

to a different cosmological evolutions, i.e. a different Wd for a given Qdc and Qcc- 

The stability of this fixed point can be investigated by linearizing eqs. (fT5|) - (ITSll 

around the fixed point and studying how the fluctuations around the fixed point evolve 

in time. If the amplitude of the fluctuations decreases in time, the fixed point is a stable 

fixed point or attractor. Linearizing eqs. ( fTSl) - ( fTSj) around {Qdc,^cc,0,0), we get 

(5^:^ = ^ - iidc - 3^ - 3Adr]dc s^d - sx^iidcSn, + ^],c5a, (21) 

5^^'c = -1 - 2-^ + SAcf^dc U(^c + 3Ad - f^cc - 3r]cc^ - SArfl^cc SQd + ^^cc^^, 

(22) 

5n'^ = ~2il + ^\ 5nr, (23) 

6ni = -il + 2^ (24) 

where 5fla = — ^ac denote fluctuations around the fixed point. The above linear 
equations can be written as 

6Q'^ = M^p6Qp, (25) 

where a and /3 run over the 4 species. The eigenvalues of the matrix M govern how 
the amplitude of the fluctuations around the fixed point changes with time. The fixed 
point is a stable, saddle or unstable point if all the eigenvalues are negative, some of 
eigenvalues are positive or all the eigenvalues are positive, respectively. (Since we are 
dealing with real eigenvalues.) The eigenvalues of the matrix M are 

qV2 q3/2 

Ai = ^ - (1 + 3Ad - 3Ac)l^dc - 3^, 
c 

A2 = A4 = -l-2^, A3 = -2 1 + ^1. (26) 
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Hence, the fixed point {Qdc, ^cc, 0, 0) is a stable point when Ai < and a saddle point 
when Ai > 0. Using eq. f|T9l) . Ai can be written as 



It is not easy to determine the sign of Ai in general. Thus, we will determine it in 
particular cases in the next section. 

3. The attractor of cosmic evolution 

We now consider the fixed point and stability of the cosmic evolution around the present 
epoch. For simplicity, we first consider the case where A^ = Ac = In this case, eq. 
(fT9|) becomes 

1 + 2^ = ^. (28) 

C i ice 

Since c > 0, eq. fl25]) implies that 6^ > This is the lower limit of 6^ for the 

existence of the fixed point. From eqs. ffT2l) and (1281) . it is easy to show that 



The equation of state parameter of the universe is defined as 

Wa^a, (30) 



Ptotal 
W = 



rtotal , ■ 

a=r,b,m,a 

where Wa = Pa/ Pa- Since Wb = Wc = and Qr can be neglected at late time, we have 
w = Wd^d and therefore. 

w = - — . (31) 

From the lower bound of 6^, we get w < 1/3. This means that the fixed point corresponds 
to cosmic acceleration. Moreover, one can see that the universe will be in the phantom 
phase, i.e., w < —1, if 6^ > ficc or equivalently O]/^^ > c. Since the observations seem to 
indicate that Wd < —1 today, we find the value of 6^ that makes Wd < —1 at the fixed 
point. It follows from eq. fl29l) that Wd < —1 if b"^ > ^cc^dc- We now check the stability 
of the fixed point. In this case, eq. fl271) becomes 

It can be seen that it is not easy to find a point at which Ai changes sign. However, 
Ai is negative for any 6^ or c if Qdc > 1/2, i.e. Qdc > ^cc- This means that the fixed 
point at which Qdc > ^cc is a stable fixed point or attractor. Since the fixed point will 
occur when Q is positive and = = 0. At the present epoch Qr is small and can be 
neglected and with positive Q the ratio flb/^c decreases with time, so the fixed point 
{fide, flee, 0, 0) will be reached in the future. From eqs. (fTSi) - (fTSjl . one can see that if 6^ 
is large, Qb can decrease quickly with time compared with Qc- As a result, the fixed point 
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can be reached quickly near the present epoch. However, the CDM dominated epoch 
will disappear and the baryon fraction will be larger than the CDM fraction in the last 
scattering epoch due to the rapid decrease of pb/pc with time. This is excluded by the 
observed peak height ratio of the CMB power spectrum. Thus, based on observations, 
the fixed point cannot be reached near the present for this case. We will see in the next 
section that small 6^ or equivalently small flee is required by observations. Substituting 
flee = 1 — f^de into eq. fl28l) . we get 

q1/2 q3/2 

1 - 36^ + - fide - = (33) 

c c 

The above third degree polynomial for fl]/^^ will have one positive real root if 1 > 3b^. 
This root is the previously considered fixed point. In contrast, if 1 < 36^ the above 
equation will have one additional real root. This root gives another fixed point at 
smaller fl^e- It is not hard to show that this fixed point is not stable and we will not 
consider it here. 

Let us now consider the case where 7^ Ac- In this case, eq. ( HM becomes 
fl^^^ 3 

1 + 2^ = — (Arf^],, + A,fi,,) . (34) 

C i ice 

Hence, the fixed point occurs when Xdfldc + Kflce > ^cc/3. Using eqs. (fT2i) and 
we get 

Arff^dc + Acf^ cc /nr\ 



and therefore 



cc / o /- \ 

^ = a • ^^^^ 



^ce 



Similarly to the case when A^ = Ac, one can show that the fixed point occurs when 
Wd < —l/{3flde) OT w < —1/3. Since the terms Xdfldc and Ac^^cc from the interaction 
dominate at different times, the parameters A^ and Ac can be chosen such that the baryon 
fraction is smaller than the CDM fraction during the last scattering epoch although flee 
need not be very small. Thus, in this case, the attractor can be reached faster than the 
case of Ac = A^. From eq. (1271) . we see that the fixed point is the stable fixed point 
when fide > f^ce and Ac is not much larger than A^. 



4. The observational constraints 



We now check the range of parameters for which the cosmic evolution has an attractor 
for compatibility with observations. We constrain the parameters of the interacting 
holographic dark energy using the latest observational data of SNIa [16] combined with 
the CMB shift parameter derived from three-year WMAP [17] and the baryon acoustic 
oscillations (BAG) from SDSS LRG [18]. 
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The SNIa observations measure the apparent magnitude m of a supernova and 
its redshift z. The apparent magnitude m is related to the distance modulus and 
luminosity distance di of the supernova by 

lj{z) = m{z) - M = 5 logio(c^L(^)/Mpc) + 25, (37) 

where M is the absolute magnitude of the supernova. For flat space time, the luminosity 
distance is given by 

dz 

d,{z) = H,\l + z) (38) 

Here, i^o is the present value of the Hubble parameter and E{z) = H{z)/Hq. To 
constrain the holographic dark energy model, we perform fit for the parameters 
Xd, Ac, c, flco, ^60) where the subscript denotes the present value. For Qro, we compute 
its value from the CMB and Neutrino temperatures. Using eq. ([2]), one can compute 
Qdo from Qco, and ^2^0- We have to include radiation in our consideration because 
it should not be neglected when we compute CMB shift parameter. For the SNIa data, 
the parameter Hq is the nuisance parameter which needs to be marginalized out. Since 
the fit of holographic dark energy to SNIa data is sensitive to Hq [19] , we have to add 
constraints from other observations to improve the fit. For this reason, we include the 
constraints from CMB shift parameter and BAO measurement in our analysis. 

The CMB shift parameter is a quantity derived from CMB data that has been 
shown to be model independent, so it can be used to constrain cosmological models. 
The CMB shift parameter is defined as [20] 



R-^Sl ^. (39) 

where Qmo = ^co + ^w and 2;cmb = 1089 is the redshift at recombination. The estimated 
value of R from 3-year WMAP data is 1.70±0.03. From the SDSS data, the measurement 
of the BAO peak in the distribution of SDSS LRG can be used to derive the model 
independent parameter A which is defined as [H] 



A = nli^zBAo) 



-1/3 



-^BAO 



dz 



z 12/3 



(40) 



.^BAO Jo ^(^). 

where zbao = 0.35. The estimated A is A = 0.469(ns/0.98)"°-^^ ± 0.017. According to 
the 3-year WMAP data, the scalar spectral index is chosen to be Us = 0.95. 

We first consider the case where Xd = Xc = b"^ ■ For simplicity, we start by neglecting 
baryons in our consideration. Hence the attractor becomes {Qd,^c,^r) = (^dc,^cc,0), 
where Qdc and Qcc satisfy eq. ( l28i) . The region of b^ and Qcc for which the fixed point 
exists and the 99.7% confidence regions from the combined constraints are shown in 
figure 1. From this figure, one sees that a small 6^ is required by observations, so that the 
attractor that satisfies the observational constraints occurs at low flee- This implies that 
the attractor cannot occurs at present at which ~ 0.3. In the future, Qc can become 
small, so that the fixed point can exist. Nevertheless, if the attractor is reached in the 
far future, the cosmic coincidence may not be alleviated because the present value of flc 
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and Wd may not satisfy the observational constraints. In the following consideration, we 
will study how the ratio r evolves during the present epoch to check the possibility of 
alleviating the cosmic coincidence. Since the physical fixed point cannot occur if 6^ < 0, 
we perform another fit by supposing that 6^ > 0. The result is shown in figure 1. It can 
be seen that the above conclusions for the case of arbitrary b"^ are also valid for the case 
of b^ > 0. 



0.08 



0.06 



0.04 




0.4 



0.4 



Figure 1. A region of parameters and fie in which the cosmic evolution has a late 
time attractor (the yellow regions above the dashed lines), and the 99.7% confidence 
levels from the combined analysis of SNIa data , CMB shift parameter and BAO 
measurement. For the yellow region, f2c refers to ilcc, but Qc refers to Oco for the 
confidence levels. The upper panel shows the case of arbitrary b^, while the lower 
panel shows the case where b^ > 0. The blue regions represent the confidence regions 
for the case that includes baryons and prior fibo = 0.047±0.006, while the thick dotted 
lines represent the confidence levels for the case where baryons are neglected. 



The situation changes a bit when we include baryons in our consideration. From 
the previous section, we know that the cosmic evolution reaches the attractor at 
{fld,^c,^b,^r) = (^dc, ^cc, 0, 0). Siuce the present value of Qb does not vanish, this 
attractor cannot occur at present. However, the attractor can occur in the future because 
the ratio Qb/^c decreases with time due to the positive Q. It can be seen that the ratio 
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^Ib/^c decreases faster when 6^ increases. According to the observational constraints, 
a small 6^ is also required in this case, so that the attractor is slowly reached in the 
future. In order to perform a fit for this case, we use a prior flbo = 0.047 ± 0.006 from 
the 1-year WMAP data [21] because QbQ cannot be constrained very well by using only 
SNIa data, CMB shift parameter and BAO measurement. It can be seen from figure 
1 that the shape of the confidence contours does not change much when we include 
non-interacting baryons in our consideration. Obviously, the contours move to the left 
when the amount of Qbo increases. This is because the amount of Qco decreases. 



0.5 




0.3 



0.14 



0.12 



0. I - 



0.08 - 
0.06 - 




Figure 2. A region of parameters Ad and Ac for which the cosmic evolution has a late 
time attractor (the yellow regions above the dashed lines), and the 99.7% confidence 
levels from the combined analysis of SNIa data , CMB shift parameter and BAO 
measurement. The upper panel shows the case of arbitrary \^ and Ac, while the lower 
panel shows the case where Ad, Ac > 0. The blue regions represent the confidence 
regions for the case that includes baryons and prior fito = 0.047± 0.006, while the thick 
dotted lines represent the confidence levels for the case where baryons are neglected. 
The dashed lines are plotted by setting ficc = 0.27, but the thick dashed line in the 
lower panel is plotted by setting ilcc — 0-05. We note that the region above the thick 
dashed line also represents the region for which the cosmic evolution has attractor. 

We now consider the case of arbitrary Ad and Ac- We also perform a fit for the 
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case with baryons and without baryons. For the case where baryons are neglected, it 
follows from figure 2 that the attractor can occur at present (n^c ~ 0.3) for a narrow 
range of A^^ and Ac- Nevertheless, a small Qcc is required by observations if we restrict 
Arf and Ac to be positive. This implies that the attractor must occur in the future. 
Similarly to the case of Xd = Xc = b^, the attractor cannot be reached at present if 
baryons are included in the consideration. Nevertheless, for suitable values of Ac and 
Xd which satisfy the observational constraints, the attractor in this case can be reached 
faster than the case of Ac = A^. 
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Figure 3. The upper panel shows the evolution of r (solid and dashed lines) and 
r/{rH) (long dashed and dotted lines), while the lower panel shows the evolution of 
Vtc (solid and dashed lines) and fid (long dashed and dotted lines). The solid and long 
dashed lines correspond to the case where the value of 6^, c and Vt^Q equal to their best 
fit value, while the dashed and dotted lines correspond to the case where 6^ is chosen 
such that Occ is close to OcO,i-e- ^cc = 0.27. 

In order to check whether the cosmic coincidence can be alleviated for this form of 
interaction, we study the evolution of r during the present epoch. Since the evolution 
of r for the case with baryons and without baryons have nearly the same feature, we 
consider only the latter case. We first consider the case where 6^ = Ac = A^. Setting 6^, 
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c and flco equal to their best fit value, i.e. 6^ = —0.004, c = 0.84 and Qco = 0.3, the 
evolution of r and r' = r/{rH) is plotted in figure 3. From figure 3, we see that \r'\ > 1 
during the present epoch because Qco is quite different from Qcc- Due to a negative 6^, 
r and Qc become negative and consequently reach the attractor at late time. Recall 
that fic > at attractor if 6^ > flcc/-i- To solve the coincidence problem, the ratio r 
should vary slowly during the present epoch such that \r'\ < 1 today The present 
value of \r'\ will decrease if the value of Qcc gets closer to the value of Qco, i-e. the 
cosmic evolution reaches the attractor near the present. According to the yellow region 
in figure 1, the value of Qcc will increase and get closer to Qco if increases. However, 
\r'\ during the present epoch will be smaller than 1 only when 6^ is larger than the 
observational bound. The evolution of r for the case where Qcc is close to Qco is shown 
in figure 3. In this case, the value of c and Qco is the best fit value, while the value 
of is larger than the observational bound but is inside the yellow region in figure 1. 
From figure 3 we see that \r'\ < 1 during the present epoch in this case but the amount 
of flc during the early time is too small. For this reason, this case does not satisfy the 
observational constraints. Hence, in the case of Ac = A^, the coincidence problem is not 
really alleviated for this form of interaction terms. 

Next, we consider the case where Ac ^ \d- We first set c, Ac, A^ and Qco equal to 
their best fit value, i.e. c = 0.37, Ac = 0.08, A^ = —0.45 and Qco = 0.27. This choice of 
the parameter value is outside the yellow region in figure 2. Since the interaction term 
Q is dominated by SHXdfld during the present epoch, Q becomes negative at present. 
Hence, r and also Qc become negative and consequently reach the attractor at late time. 
Of course, the attractor in this case does not correspond to the physical attractor. The 
evolution of r and Qc is shown in figure 4. Similarly to the case of Ac = A^, the present 
value of \r'\ for this choice of parameters value is larger than 1 because Qcc is quite 
different from Q^o- We now keep Qco fixed and choose the new value of parameters c, 
Ac and Xd such that it satisfies the observational constraints and the cosmic attractor 
occurs at Qcc = 0.24 near the present epoch. This choice of the parameters value is 
inside the intersection between the yellow region and the confidence region in figure 2. 
From figure 4, we see that the present value of |r'| is smaller than 1 and the attractor 
is reached near the present epoch for this choice of the parameters value. Based on soft 
coincidence idea, the coincidence problem can be alleviated in this case. The coincidence 
problem is better alleviated if Qcc = ^co or equivalently if the cosmic attractor is reached 
at the present. Unfortunately, the parameter values that make flee = ^co do not satisfy 
the observational constraints. The evolution of r and Qc for this case is shown in figure 
4. From the figure, it is clear that this case is ruled out by observations. We note that 
the present value of \r'\ will increase if we include baryons in the above consideration. 

Finally, we estimate whether \r'\ can be smaller than 1 at present if we constrain the 
parameters of dark energy by SDSS matter power spectrum. Here, we will not perform 
a complete fit for this dataset because the nature of the perturbations in holographic 
dark energy is not yet completely understood. To write down the evolution equations for 
the density perturbation and compute the matter power spectrum for this interacting 
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Figure 4. The upper panel shows the evolution of r (solid, dashed and long dashed 
dotted lines) and f/{rH) (long dashed, dotted and dashed dotted lines), while the 
lower panel shows the evolution of fie (solid, dashed and long dashed dotted lines) and 
fid (long dashed, dotted and dashed dotted lines). The solid and long dashed lines 
correspond to the case where the value of Ac, A^, c and f2co equal to their best fit value, 
the dashed and dotted lines correspond to the case where c, Ac and A^ are chosen such 
that Vice is close to f2co,i-e. Vice = 0.24, and the long dashed dotted and dashed dotted 
lines correspond to the case where c, Ac and A^ are chosen such that Vice = Vlco- 



dark energy model, we use the assumption below. From the results in [22], we suppose 
that the perturbation in holographic dark energy can be neglected when kRe/a ^ 1, 
where k is the wavenumber of the perturbation modes. Using eqs. ([1]) and ([6]) and 
supposing that radiation can be neglected during matter domination, one can show 
that H~^/Re = y/Qd/c, i.e. the Hubble radius is smaller than the event horizon if 
\/Tld < c. Since we use the data from SDSS which measures the matter power spectrum 
on scales smaller than the Hubble radius, we neglect the perturbation in holographic 
dark energy in the calculation of the matter power spectrum. We write the perturbed 
interaction term using the formulas in [231 El] , so that the evolution equations for the 
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perturbation in CDM are 
dAr 



drj 
dV, 
drj 



- kV^ + Ae + A 



Pd 



drj 



Ac + A, 



pd 



pd 



= --HV^ + k-^, 



(41) 



where Ac and Vc are the gauge-invariant density contrast and velocity perturbation of 
CDM, \1/ and $ are the metric perturbation, Ti = a~^{da/dri) and rj is the conformal 
time. We solve the above equations, compute the matter power spectrum and compare 
the obtained matter power spectrum with SDSS data using CMBEASY [25]. By 
checking the value of x^i '^6 have found that the best fit parameters for this dataset 
are different from the best fit parameters from the observational constraints on the 
homogeneous universe. Instead of searching for the best fit parameters for this dataset, 
we roughly check the viability of the parameters by comparing the matter power 
spectrum of the considered models with the matter power spectrum of ACDM model 
whose parameters are taken from the best fit value of 3-year WMAP and SDSS data 
|26] . In figure 5, we plot the fractional difference of the matter power spectrum between 
interacting holographic dark energy and ACDM model. From this figure, we see that for 
suitable ranges of dark energy parameters, \r'\ can be smaller than 1 at present and the 
difference of the matter power spectrum can be smaller than the error for the matter 
power spectrum of SDSS data. This implies that the alleviation of the coincidence 
problem by this interacting dark energy model is not excluded by SDSS data. 




0.001 



0.01 0.1 

fc//i(Mpc-i) 



Figure 5. The fractional difference of the matter power spectrum \5P{k) / P{k)\ 
between interacting holographic dark energy and ACDM model. For this plot, the 
parameters of dark energy are chosen such that \r'\ k, 0.9 at present. The fractional 
error for the matter power spectrum \6P{k)/ P{k)\ = error / spectrum of SDSS data 
is represented by dots. 
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5. Conclusions 

For the interacting holographic dark energy model, we study the fixed points and their 
stabiliby, and compare a range of model parameters for which attractor exists with the 
99.7% confidence levels from the combined analysis of SNIa data, CMB shift parameter 
and BAO measurement. Neglecting baryons, the observational constraints require that 
the value of at the attractor point must be small if = Ac or A^, Ac > 0. This 
implies that the cosmic evolution will reach the attractor point in the future when flc 
becomes small. In this case, r cannot be slowly varying during the present epoch and the 
cosmic attractor cannot be reached near the present. Hence, the coincidence problem 
is not really alleviated for this case. However, if A^ and Ac are allowed to be negative, 
the cosmic evolution can reach the attractor near the present epoch for a narrow range 
of Xd and Ac. Therefore, the coincidence problem is possible to alleviate in this case. 
Including baryons in our consideration, the attractor of the cosmic evolution cannot 
occur at present due to the non- vanishing baryon fraction. According to observations, 
the fixed point in this case is possible only when Q is small and positive. Hence, the fixed 
point will be slowly reached in the future. These results indicate that for the interacting 
holographic dark energy model with the interaction terms considered here, the cosmic 
coincidence problem cannot be alleviated very well. We also briefly considered the 
constraint from SDSS matter power spectrum on the dark energy parameters. We have 
found that the parameters ranges that lead to the alleviation of cosmic coincidence are 
allowed by SDSS data. 
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